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[i] Ozone (0 3 ) precursor emissions influence regional and global climate and air quality 
through changes in tropospheric 0 3 and oxidants, which also influence methane (CH 4 ) 
and sulfate aerosols (S0 4 ). We examine changes in the tropospheric composition of 0 3 , 
CH 4 , SO 4 and global net radiative forcing (RF) for 20% reductions in global CH 4 burden 
and in anthropogenic 0 3 precursor emissions (NO x , NMVOC, and CO) from four 
regions (East Asia, Europe and Northern Africa, North America, and South Asia) using 
the Task Force on Hemispheric Transport of Air Pollution Source-Receptor global 
chemical transport model (CTM) simulations, assessing uncertainty (mean ± 1 standard 
deviation) across multiple CTMs. We evaluate steady state 0 3 responses, including 
long-term feedbacks via CH 4 . With a radiative transfer model that includes greenhouse 
gases and the aerosol direct effect, we find that regional NO x reductions produce 
global, annually averaged positive net RFs (0.2 ± 0.6 to 1.7 ± 2 mWm~ 2 /Tg N yr -1 ), 
with some variation among models. Negative net RFs result from reductions in 
global CH 4 (—162.6 ± 2 mWm -2 for a change from 1760 to 1408 ppbv CH 4 ) and 
regional NMVOC (—0.4 ± 0.2 to —0.7 ± 0.2 mWm 2 /Tg C yr ') and CO emissions 
(—0.13 ± 0.02 to —0.15 ± 0.02 mWm -2 /Tg CO yr -1 ). Including the effect of 0 3 
on C0 2 uptake by vegetation likely makes these net RFs more negative by —1.9 to 
—5.2 mWm -2 /Tg N yr -1 , —0.2 to —0.7 mWm -2 /Tg C yr -1 , and —0.02 to —0.05 mWm -2 / 
Tg CO yr . Net RF impacts reflect the distribution of concentration changes, where 
RF is affected locally by changes in S0 4 ~, regionally to hemispherically by 0 3 , and 
globally by CH 4 . Global annual average S0 4 responses to oxidant changes range from 
0.4 ± 2.6 to — 1.9 ± 1.3 Gg for NO x reductions, 0.1 ± 1.2 to —0.9 ± 0.8 Gg for NMVOC 
reductions, and —0.09 ± 0.5 to —0.9 ± 0.8 Gg for CO reductions, suggesting additional 
research is needed. The 100 -year global warming potentials (GWP 100 ) are calculated 
for the global CH 4 reduction (20.9 ±3.7 without stratospheric 0 3 or water vapor, 

24.2 ± 4.2 including those components), and for the regional NO x , NMVOC, and CO 
reductions (-18.7 ± 25.9 to -1.9 ± 8.7 for NO x , 4.8 ± 1.7 to 8.3 ± 1.9 for NMVOC, 
and 1.5 ± 0.4 to 1.7 ± 0.5 for CO). Variation in GWP 10 o for NO x , NMVOC, and CO 
suggests that regionally specific GWPs may be necessary and could support the inclusion 
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of O3 precursors in future policies that address air quality and climate change simultaneously. 

Both global net RF and GWP100 are more sensitive to NO x and NMVOC reductions 
from South Asia than the other three regions. 

Citation: Fry, M. M., et al. (2012), The influence of ozone precursor emissions from four world regions on tropospheric 
composition and radiative climate forcing, J. Geophys. Res., 117 , D07306, doi: 10. 1029/201 1JD0 17 134. 


1. Introduction 

[ 2 ] Tropospheric ozone (0 3 ), methane (CH 4 ), and aerosols 
make important contributions to the global mean radiative 
forcing (RF) of climate [ Forster et al., 2007; Ramaswamy 
et al., 2001], Here we aim to quantify the net RF of these 
species due to regional changes in 0 3 precursor emissions, 
across an ensemble of global chemical transport models 
(CTMs). We define net RF as the net (incoming minus 
outgoing) change in irradiance (solar and infrared) at the 
tropopause between a base and perturbed state (in Watts per 
square meter [W nr -2 ]) after allowing stratospheric tem- 
peratures to readjust [Forster et al, 2007]. The contribution 
of changes in tropospheric O3 to the global mean RF since 
preindustrial times is an estimated 0.35 ± 0.15 Wm~ 2 , 
which is approximately 21% of the RF due to changes in 
carbon dioxide (C0 2 ) [Forster et al., 2007]. Changes in CH 4 
have contributed approximately 0.48 Wm" 2 , while those in 
sulfate aerosols (S0 4 - ) have contributed approximately 
—0.4 ± 0.2 Wm~ 2 (direct effect only) [Forster et al., 2007]. 

[3] Changes in O3 precursor emissions (nitrogen oxides 
[NO x ], non-methane volatile organic compounds [NMVOC], 
carbon monoxide [CO], and CH 4 ) affect the abundance of 
gaseous species (0 3 and CH 4 ), and aerosols via changes in 
the availability of atmospheric oxidants (hydroxyl radical 
[OH], hydrogen peroxide [H 2 0 2 ], 0 3 ) [Pham et al., 1995; 
Unger et al., 2006; Shindell et al., 2009; Leibensperger 
et al., 2011], These perturbations in turn influence the RF 
due to 0 3 and CH 4 and inorganic aerosol-phase species 
[Ming et al, 2005; Unger et al., 2006; Naik et al, 2007; 
Shindell et al., 2009]. O3 precursors also affect organic 
aerosols, including the formation of secondary organic 
aerosols (SOA) [Carlton et al., 2010], but the resulting RF 
remains to be quantified. O3 decreases plant growth and 
hence reduces the removal of C0 2 from the atmosphere 
[Felzer et al., 2007; Sitch et al, 2007; Collins et al., 2010], 
while NO x emissions influence nitrogen deposition and the 
subsequent uptake of C0 2 in terrestrial and oceanic eco- 
systems [Holland and Lamarque, 1997; Duce et al., 2008]. 
Because of these influences, actions to control O3 precursor 
emissions affect both air quality and global climate [Fiore 
et al., 2002; Dentener et al., 2005; West et al., 2007]. 

[ 4 ] Past studies have shown that both regional and global 
reductions in NO x surface emissions produce an overall 
positive RF from global CH 4 increases via decreases in OH, 
which outweigh the negative forcing from tropospheric O3 
decreases. The magnitude of forcing, however, depends on 
the location or sector of emission reduction [Fuglestvedt 
et al, 1999; Wild et al, 2001; Berntsen et al., 2005; Naik 
et al., 2005; West et al., 2007; Derwent et al., 2008], 
In contrast, CO, NMVOC, and CH 4 reductions contribute 
an overall negative RF by decreasing tropospheric O3 and 
increasing OH, leading to global CH 4 decreases [Prather, 
1996; Wild et al., 2001; Fiore et al., 2002; Naik et al., 
2005], Global anthropogenic CH 4 emission reductions were 


shown to produce the most negative RF of the O3 precursors, 
mainly due to direct reductions in CH 4 forcing [Fiore et al., 
2002; Shindell et al., 2005; West et al, 2007]. 

[5] Because of the short atmospheric lifetimes of O3, 
aerosols, and their precursors (apart from CH 4 ), studies of 
regional O3 precursor reductions show that air quality and 
RF effects depend strongly on the geographical location 
of emissions [Fuglestvedt et al., 1999; Berntsen et al., 2005; 
Naik et al., 2005]. The dependence on location, however, 
has made it difficult to include O3 precursors in emissions 
trading schemes or international climate agreements [Rypdal 
et al., 2005], and to evaluate the co-benefits of actions to 
reduce 0 3 for slowing global climate change. Studies have 
shown substantial dependence of O3 concentrations and RF 
on the region of NO x emissions [Fuglestvedt et al., 1999; 
Berntsen et al., 2005; Naik et al., 2005 ; Derwent et al., 2008; 
Fuglestvedt et al., 2010], where O3 and RF are more sensi- 
tive to NO x reductions from tropical regions, yet positive 
global annual average net RFs result from NO x reductions in 
each of nine world regions [Naik et al., 2005], While 
NMVOCs typically have atmospheric lifetimes comparable 
to NO x , CO has a longer atmospheric lifetime of 1-3 months 
[Seinfeld and Pandis, 2006], suggesting that the influence of 
CO on 0 3 and OH may be less dependent on reduction 
region [Berntsen et al., 2005; Rypdal et al., 2005]. Apart 
from Naik et al. [2005], who also examined combined 
reductions in NO x , CO, and NMVOCs from three regions, 
and Berntsen et al. [2005], who evaluated CO changes from 
two regions, less attention has been placed on regional CO 
and NMVOC emissions and their effects on tropospheric O3, 
CH 4 , and S0 4 ~ and global and regional net RF. In response 
to nonuniform forcings, some regional climate responses are 
also sensitive to the location and distribution of O3, aerosols, 
and RF [Shindell and Faluvegi, 2009], 

[6] Here we investigate the effects of a 20% reduction in 
global CH 4 abundance and 20% reductions in anthropogenic 
emissions of NO x , NMVOC, and CO, individually and 
combined, from four world regions on tropospheric 0 3 , 
CH 4 , and S0 4 - concentrations and on the resulting distri- 
bution and magnitude of global net RF for all precursor- 
region pairs. We use results from an ensemble of global 
CTMs that participated in the Task Force on Hemispheric 
Transport of Air Pollution (TF HTAP) multimodel inter- 
comparison study of Source-Receptor (SR) sensitivity [Fiore 
et al., 2009], which allows for an assessment of uncertainty 
as the spread across CTMs. We examine the regional depen- 
dency of RF and global warming potential (GWP) to indi- 
vidual precursors by comparing estimates across the four 
regions of reduction. 

2. Methodology 

[ 7 ] We use the results from an ensemble of 1 1 global 
CTMs [Task Force on Hemispheric Transport of Air 
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Table 1. HTAP Source-Receptor Sensitivity Simulations, 
Where the Four Regions of Reduction are East Asia, Eur- 
ope, North America, and South Asia for SR3 Through SR 6 


Scenario Description 


SRI 

Base simulation 

SR2 

20% reduction in global CH4 
mixing ratio 

SR3 

20% reduction in regional 
NO x emissions 

SR4 

20% reduction in regional 
NMVOC emissions 

SR5 

20% reduction in regional 
CO emissions 

SR6 

20% reduction in regional NO x , 

NMVOC, CO, and aerosol emissions 


Pollution, 2010] to evaluate changes in the tropospheric 
distributions of 0 3 , CH 4 , and S0 4 - for each reduction sce- 
nario. The net RF is calculated from the simulated changes 
in 0 3 , CH 4 , and S0 4 - using the standalone radiative transfer 
model (RTM) developed by the National Oceanographic and 
Atmospheric Administration (NOAA) Geophysical Fluid 
Dynamics Laboratory (GFDL) [ Schwarzkopf and Ramaswamy, 
1999; GFDL Global Atmospheric Model Development Team 
(1 GFDL GAMDT), 2004], We estimate the RF due to S0 4 “, 
as S0 4 _ responds to oxidant changes, considering only its direct 
effects on radiation [Naik et al., 2007]; few CTMs reported 
other aerosol species for all reduction scenarios. Although we 
account for variability across die CTM ensemble in the O 3 , 
S0 4 ~, and CFI 4 estimates, we likely underestimate the uncer- 
tainty in net RF by using only a single RTM. 

2.1. HTAP CTM Simulations 

[8] The SR simulations performed by each CTM are out- 
lined in Table 1. We analyze O 3 and CH 4 results from 11 
CTMs and S0 4 _ results from four CTMs (Table 2). Each 
CTM utilized its own emissions inventory and prescribed 
meteorological fields for the year 2001 [Fiore et al., 2009]. 
Anthropogenic emissions of NO x , NMVOCs, CO, and all 
precursors combined were reduced by 20% in each of four 
world regions: East Asia (EA), Europe and Northern Africa 
(EU), North America (NA), and South Asia (SA). For CH 4 , 
the present-day abundance (1760 parts per billion by volume 
[ppbv]) was imposed in all simulations except for the CH 4 


Table 2. Global CTMs Used for Multimodel Mean 0 3 , CH 4 , 
and S0 4 ~ Estimates 


Model 

Institution 

CAMCHEM-331 lml3 a 

NCAR, USA 

FRSGCUCI-vOl 

Lancaster University, UK 

GISS-PUCCINI-modelE 

NASA GISS, USA 

GMI-v02f 

NASA GSFC, USA 

INCA-vSSz 

LSCE, France 

LLNL-IMP ACT -T5 a a 

LLNL, USA 

MOZARTGFDL-v2 a 

GFDL, USA 

MOZECH-V16 

FZ Juelich, Germany 

STOC-HadAM3-v01 

University of Edinburgh, UK 

TM5-JRC-cy2-ipcc-vl a 

JRC, Italy 

UM-CAM-vOl 

University of Cambridge, 
UK/NIWA, NZ 


a The four global CTMs used for multimodel mean SO4 estimates. 


control simulation (SR2), where global CH 4 abundance 
was decreased by 20% to 1408 ppbv. All simulations were 
performed for a full year (2001), after an initial spin-up of 
at least six months [Fiore et al., 2009], 

[ 9 ] The multimodel mean ± 1 standard deviation changes 
in the anthropogenic emissions of NO x , NMVOCs, and CO, 
across 1 1 CTMs, are presented in Table 3. There is consid- 
erable variability in the emission reduction magnitudes 
across CTMs. Coefficients of variation (CVs) (standard 
deviation/mean) are lowest for NO x emissions from EU, 
NA, and SA, while there is more variability in NMVOC and 
CO emissions from the four regions, consistent with the 
comparison of current global emission inventories by 
Granier et al. [2011]. 

[ 10 ] Previous publications based on the HTAP SR 
experiments have emphasized the effects of long-range 
transport on surface O 3 , and other components, and com- 
parisons with observations [Sanderson et al., 2008; Shindell 
et al., 2008; Fiore et al, 2009; Jonson et al., 2010; 
Reidmiller et al., 2009]. Fiore et al. [2009] found that the 
HTAP ensemble mean surface O 3 concentrations compared 
well with observations over EU for the year 2001, but 
overestimated measurements by more than 1 0 ppb during the 
summer and fall over the eastern United States and Japan. 
Jonson et al. [2010] compared simulated vertical O 3 profiles 
with observed ozonesonde profiles, finding that the spread in 
CTM results (and their over and underestimation of 0 3 
soundings) increases in the spring and summer with more 
active chemistry. In the winter and spring, seasonal averages 
for most CTMs were within 20% of sonde measurements in 
the upper and middle troposphere. Simulated S0 4 _ con- 
centrations at the surface for the base simulation (SRI) also 
have been compared to observations (M. Schulz, personal 
communication, 2011, preliminary results available at http:// 
aerocom.met.no/cgi-bin/aerocom/surfobs_annualrs.pl), where 
the results show that the CTMs are generally realistic. 

[ 11 ] Short-lived 0 3 precursors affect tropospheric 0 3 
within hours to weeks after their emission; however, they 
also affect OH, which influences the lifetime of CH 4 and in 
turn, O 3 in the long term [Prather, 1996; Wild et al., 2001; 
Berntsen et al, 2005; Naik et al, 2005]. Global CH 4 chan- 
ges were calculated by Fiore et al. [2009], based on the CH 4 
loss by tropospheric OH diagnostic reported for each CTM 
and SR3 through SR6, relative to the fixed CH 4 abundance 
of 1760 ppbv. Long-term 0 3 responses are then calculated in 
each grid cell, following West et al. [2007, 2009b], by 
scaling the change in O 3 from the CH 4 control simulation 
(SR2 minus SRI) to the calculated global CH 4 change for 

Table 3. Multimodel Mean ± 1 Standard Deviation Reductions in 
the Anthropogenic Emissions of NO x , NMVOC, and CO (20% 
of Total Anthropogenic Emissions) Among the 1 1 HTAP CTMs 
Used Here 3 


Region 

NO x (Tg N a- 1 ) 

NMVOC (Tg Ca^ 1 ) 

CO (Tg a- 1 ) 

EA 

EU 

NA 

SA 

1.17 ± 0.24 (0.20) 
1.48 ± 0.14 (0.09) 
1.48 ± 0.10 (0.07) 
0.46 ± 0.04 (0.09) 

3.13 ± 1.24 (0.40) 
3.77 ± 1.88 (0.50) 
3.11 ± 1.34 (0.43) 
1.94 ± 0.63 (0.33) 

25.58 ± 7.25 (0.28) 
15.40 ± 3.26 (0.21) 
19.69 ± 3.66 (0.19) 
15.82 ± 3.74 (0.24) 


“Coefficient of variation (CV = standard deviation/mean) is in 
parentheses. Emissions for the individual CTMs are provided in Tables A2 
and A3 of Fiore et al. [2009], 
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Figure 1. Global annual average changes in full (blue) and upper (yellow) tropospheric 0 3 burden (Tg) 
at steady state (perturbation minus base), where the upper troposphere is from 500 hPa to the tropopause, 
for the HTAP ensemble of 1 1 models, showing the median (black bars), mean (red points), mean ± 1 SD 
(boxes), and max and min (whiskers), for each precursor reduction scenario (—20% global CH 4 burden, 
and —20% regional emissions of NO x , NMVOC, CO, and combined from East Asia [EA], Europe and 
Northern Africa [EU], North America [NA], and South Asia [SA]). 


each SR simulation and CTM. We then add the long-term 0 3 
responses to the simulated short-term 0 3 responses to give 
0 3 concentrations at steady state. 

[ 12 ] Since the HTAP CTMs were typically not designed 
to model stratospheric chemistry, we use the same three- 
dimensional monthly mean 0 3 concentrations in the 
stratosphere and merge these with calculated steady state 
(short-term + long-term) tropospheric 0 3 concentrations for 
each simulation and CTM. Stratospheric 0 3 is taken for the 
year 2001 from the AC&C/SPARC 0 3 database prepared 
for CMIP5 (Available: http://pcmdi-cmip.llnl.gov/cmip5/ 
forcing.html). Sovde et al. [2011] found that around 15% of 
the RF from 0 3 precursors is due to 0 3 changes in the lower 
stratosphere, using a single model with both standard and 
updated chemistry. Since we ignore lower stratospheric 0 3 
changes, our RF estimates may underestimate the full effect 
of 0 3 precursors. After each CTM’s 0 3 and S0 4 - results are 
interpolated to a common resolution (longitude x latitude x 
level) as required by the RTM (72 x 37 x 33 for 0 3 ; 96 x 
80 x 14 for SOI - ), the HTAP ensemble mean ± 1 standard 
deviation 0 3 and SO| - distributions are calculated in each 
grid cell and month in three dimensions, in addition to the 
ensemble mean ± 1 standard deviation global CH 4 abun- 
dances (derived from the CH 4 loss by tropospheric OH 
diagnostics). Global 0 3 , CH 4 , and S0 4 - changes are calcu- 
lated for each CTM as perturbation (SR2 to SR6) minus base 
(SRI) values. 

2.2. GFDL Radiative Transfer Model 

[ 13 ] We employ the GFDL RTM to estimate the net RF at 
steady state due to the changes in atmospheric gases (0 3 and 
CH 4 ) alone and due to combined changes in 0 3 , CH 4 , and 
S0 4 - aerosols. The GFDL RTM is a module of the GFDL 
coupled atmospheric-ocean model (AM2) and simulates 
solar and infrared radiative transfer [GFDL GAMDT, 2004; 


Naik et al., 2005, 2007], It has been used in studies of long- 
lived greenhouse gases [Schwarzkopf and Ramaswamy, 
1999] and short-lived forcing agents such as 0 3 and aero- 
sols [Naik et al., 2005, 2007; West et al, 2007; Fiore et al, 
2008; Saikawa et al., 2009]. Here the RTM is employed as 
for Naik et al. [2007] and Saikawa et al. [2009], at 144 x 
90 x 24 levels, except for the following changes. We update 
well-mixed greenhouse gas concentrations based on obser- 
vations for the year 2001 included as part of the historical 
period (1750-2005) of the CMIP5 Representative Con- 
centration Pathways (RCP) database [Meinshausen et al., 
2011] (Available: http://www.iiasa.ac.at/web-apps/tnt/RcpDb/ 
dsd?Action=htmlpage&page=download). We also update 
the solar data used by the RTM to the CMIP5 solar forc- 
ing data (Available: http://www.geo.fu-berlin.de/en/met/ag/ 
strat/forschung/SOLARlS/Input_data/CMIP5_solar_irradiance. 
html). The RTM simulations do not include the indirect 
effects of aerosols on clouds or the internal mixing of aero- 
sols. RF contributions from changes in nitrate aerosols, 
changes in stratospheric 0 3 and water vapor, changes to the 
carbon cycle via 0 3 and nitrogen deposition, and changes to 
CO 2 from CH 4 , CO, and NMVOC oxidation are also omitted 
in the RTM simulations. 

[ 14 ] The multimodel monthly mean ± 1 standard devia- 
tion 0 3 , CH 4 , and S0 4 - concentrations are used as input in 
the RTM simulations, along with meteorological fields from 
GFDL’s atmosphere model (AM2) and land model (LM2), 
sampled one day per month at midmonth for the year 2001 
to represent monthly mean conditions [Naik et al., 2005], 
Substantial variability in the S0 4 - estimates across only four 
CTMs precluded evaluating the mean — 1 standard deviation 
for S0 4 - (e.g., several grid cells had standard deviations 
exceeding the mean). We simulate the monthly mean net 
radiation fluxes for the base and perturbed cases and calcu- 
late the net RF as the difference between the perturbed and 
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Figure 2. Global annual average changes in full (blue) and upper (yellow) tropospheric 0 3 burden per 
change in emissions (Tg 0 3 /Tg emissions per year) at steady state for the individual 11 models, where 
the units of emissions are Tg N (for NO x ), Tg C (for NMVOCs), and Tg CO (for CO), showing the 
median (black bars), mean (red points), mean ± 1 SD (boxes), and max and min (whiskers) across the 
HTAP ensemble. 


base net fluxes (net shortwave minus net longwave) at the 
tropopause, after allowing stratospheric temperatures to 
readjust to radiative equilibrium [Naik et al., 2007; Saikawa 
et al, 2009]. 

3. Tropospheric Composition Changes 
3.1. Tropospheric Ozone Changes 

[is] Figure 1 shows the changes in global annually aver- 
aged steady state tropospheric 0 3 burden and its variability 
across 1 1 HTAP CTMs. Full troposphere and upper tropo- 
sphere (UT) 0 3 burdens are distinguished because 0 3 in the 
UT has a higher RF efficiency on a per molecule basis [Lacis 
et al., 1990; Wang et al, 1993; Forster and Shine, 1997]. 
For each CTM’s regridded 0 3 distributions that have been 
blended with CM1P5 stratospheric 0 3 values (section 2.1), 
the UT is defined from 500 hPa to the tropopause, where the 
tropopause is identified at the 150 ppbv 0 3 level. 

[ 16 ] The largest changes in full troposphere and UT 0 3 
burden are found for the 20% CH 4 reduction, followed 
by the 20% combined precursor reductions from NA and 
EA, respectively. However, there is considerable diversity 
among the 11 CTMs’ estimates of full troposphere and UT 
0 3 burden changes. In these 17 SR simulations relative 
to the base case, the change in UT 0 3 burden per change in 
full troposphere 0 3 burden (UT 0 3 /fiill troposphere 0 3 ) is 
largest for reductions in global CH 4 (0.36 to 0.47) and regional 
CO emissions (0.19 to 0.53), and smallest for regional 
NMVOC reductions (0.16 to 0.42) across the 11 CTMs, 
reflecting the longer lifetimes of CO and CH 4 in comparison to 
NMVOCs and NO x . UT 0 3 /full troposphere 0 3 is also largest 
for individual 0 3 precursor (NO x , NMVOC, and CO) reduc- 
tions from SA in comparison to the other regions. 

[ 17 ] To evaluate regional sensitivities, we consider how 
changes in full troposphere and UT 0 3 burdens per unit 
change in emissions vary by region for each precursor 
(Figure 2). The CTMs mostly agree that the SA NO x 
reduction produces the greatest change in full troposphere 
and UT 0 3 burden per change in emissions out of the four 


regions, which can be attributed to more rapid vertical 
mixing, stronger photochemistry, and greater sensitivity of 
0 3 to precursor emissions from the tropics and southern 
hemisphere (SH) \Kunhikrishnan et al, 2004; West et al, 
2009a]. We find less variability across the four regions in 
reducing full troposphere and UT 0 3 burden per change in 
NMVOC and CO emissions (Figure 2), but six CTMs show 
that SA NMVOC and CO reductions produce the largest 
reductions in UT 0 3 burden per change in emissions. 

[is] For each regional reduction, the greatest changes in 
steady state tropospheric total column 0 3 occur over the 
reduction region, in each individual CTM and across the 
HTAP ensemble (Figure 3). For NO x reductions from each 
region, slight increases in 0 3 burden occur in the SH, in 
contrast to overall decreases in the northern hemisphere 
(NH), due to long-term 0 3 increases via CH 4 that are glob- 
ally distributed (according to the pattern of 0 3 response to 
CH 4 ) [West et al, 2007, 2009b] (see section 3.2). As SA is 
further south than the other regions, the largest 0 3 decreases 
occur near the tropics and do not extend as far north. 
The slightly higher increases in 0 3 in the SH from the NA 
NO x reduction correspond to this region producing greater 
increases in global CH 4 (Figure 4). For NMVOC and CO 
reductions from the four regions, we find decreases in total 
column 0 3 in both the SH and NH, as sustained decreases 
in these precursors cause both short- and long-term global 
0 3 decreases. NMVOC reductions also reduce global annual 
average peroxy acetyl nitrate (PAN) burdens from the four 
regions by 1.4 ± 0.6% (SR4EA), 2.0 ± 0.8% (SR4EU), 
1.5 ± 0.5% (SR4NA), and 0.5 ± 0.2% (SR4SA), relative to 
the base, through which NMVOCs can influence the nitro- 
gen cycle and therefore long-range 0 3 . Across the 1 1 CTMs, 
we find that tropospheric PAN decreases are correlated 
to the changes in NMVOC emissions for the EA, EU, 
and NA reductions. 

3.2. Tropospheric Methane Changes 

[ 19 ] Although global CH 4 was held constant by each CTM 
in all perturbations, we analyze the changes in global 
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Figure 3. Annual average steady state tropospheric total column O 3 burden changes ( 1 0 — 2 DU) for the 
multimodel mean of 1 1 HTAP models, for each of the precursor reduction scenarios (—20% CH 4 burden, 
and —20% regional emissions of NO x , NMVOC, CO, and combined). The 4 regions of reduction 
(NA, EU, SA, EA) are outlined in red in the —20% CH 4 plot. 


tropospheric CH 4 calculated off-line for each perturbation. 
NO x reductions from all four regions increase global CH 4 
abundance via decreases in OH, while NMVOC and CO 
reductions from all four regions decrease global CH 4 
(Figure 4). These changes in CH 4 drive the long-term O3 
changes discussed in the previous section. For particular 
precursors, reductions from certain regions, e.g., CO reduc- 
tions from EA, are slightly more effective at decreasing 
global CH 4 than other regions. However, there is notice- 
able variability among the CTMs’ changes in global CH 4 
(Figure 4), which is partly explained by variance in CTM 
emissions for CO, but not for NO x and NMVOCs. 

[ 20 ] For several CTMs, some emissions perturbations had 
minimal impact on global OH, resulting in calculated steady 
state CH 4 changes of zero. In addition, the CVs of CH 4 
change are lowest in magnitude for NO x reductions (0.22 to 


0.39) and highest for NMVOC reductions (—0.40 to —1.12), 
perhaps reflecting differing NMVOC speciation and chemi- 
cal schemes among the CTMs [ Collins et al., 2002]. 

[ 21 ] We compare the change in global CH 4 from the 
combined precursor reductions to the sum of global CH 4 
changes from NO x , NMVOC, and CO reductions to assess 
the additivity of individual precursor reductions on global 
CH 4 . Based on the HTAP multimodel mean results, global 
CH 4 changes from SA and NA combined precursor reduc- 
tions are approximately 59% and 75%, respectively, of the 
sum of CH 4 changes from individual precursor reductions, 
while global CH 4 changes from EA and EU combined pre- 
cursor reductions are slightly negative in contrast to small 
positive global CH 4 changes from the sum. The three (of 11) 
CTMs that did not include reductions in SO 2 and aerosols in 
the SR 6 experiments [Fiore et al., 2009] show global CH 4 
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-20% NOx -20% NMVOC -20% CO -20% Combined 

Figure 4. Global annual multimodel changes (perturbation minus 1760 ppbv) in tropospheric CH 4 
(ppbv) for —20% regional emissions of NO x , NMVOC, CO, and combined: median (black bars), mean 
(red points), mean ± 1 SD (boxes), and max and min (whiskers) for the HTAP ensemble of 1 1 models, 
estimated directly from the CH 4 loss by tropospheric OH archived by each HTAP CTM [Fiore et al., 
2009]. Tropospheric CH 4 changes were not available from INCA-vSSz for SA 20% NMVOC reduction 
and from LLNL-IMPACT-T5a for 20% NO x reductions (EA, EU, NA, SA); these models are excluded 
from the multimodel CH 4 changes for these perturbations. 


changes from the combined precursor reductions close to 
the sum of CH 4 changes from individual NO x , NMVOC, 
and CO reductions. This suggests that deviations from 
additivity may be due to reductions in SO 2 and aerosols 
(in SR6) affecting CH 4 . 

3.3. Tropospheric Sulfate Changes 

[ 22 ] There is considerable variability and disagreement in 
the sign of S0 4 _ responses among the four CTMs evaluated 
(Figure 5). The greatest variability in global S0 4 _ burden 
across the CTMs occurs for the CH 4 reduction and for NO x 
reductions from all four regions. There is less variability 
across CTMs for NMVOC and CO reductions, but still dif- 
ferences in the sign of change. 


[ 23 ] The four CTMs analyzed here account for S0 4 - for- 
mation through two main oxidation pathways: 1) gas-phase 
oxidation of sulfur dioxide (S0 2 ) by OH, and 2) aqueous- 
phase oxidation of S0 2 by H 2 0 2 or O3 [Houweling et al., 
1998; Jeuken et al., 1999, 2001; Barth et ah, 2000; Rasch 
et al., 2000; Horowitz et al., 2003; Rotman et al., 2004; 
Tie et al., 2005]. Since CH 4 , NMVOC, and CO reductions 
increase OH concentrations, S0 2 oxidation via pathway 
(1) is expected to increase tropospheric S0 4 ~ formation. 
At the same time, increases in H 2 0 2 (occurring with 
increases in OH) are expected to increase S0 2 oxidation via 
pathway (2), but decreases in 0 3 may decrease oxidation 
(also pathway (2)). NO x reductions not only decrease O3, but 
can also decrease OH and H 2 0 2 , which leads to decreases in 



EA EU NA SA EA EU NA SA EA EU NA SA 



-20% NO x -20% NMVOC -20% CO -20% CH 4 

Figure 5. Global annual multimodel changes (perturbation minus base) in short-term tropospheric S0 4 
(Gg) for —20% CH 4 burden and —20% regional emissions of NO x , NMVOC, and CO: mean (red bars) 
and mean ± 1 SD (boxes) across the HTAP ensemble of four models. The individual model results are 
shown in black (pluses). 
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Figure 6. Annual average tropospheric total column SO4 - burden changes (jig m 2 ) for the multimodel 
mean of four HTAP models for —20% CH 4 burden and —20% regional emissions of NO x , NMVOC, 
and CO scenarios. 


SOI - formation by pathways (1) and (2). Figure 5 shows 
that the sign of SO 4 - response is not consistent across all 
four CTMs, suggesting uncertainty in the modeled effects 
of O3 precursors on oxidants, the relative importance of 
these oxidation pathways, and the lifetime and removal of 
S0 2 and SO| - . 

[24] The global distributions of tropospheric total column 
S0 4 - changes (Figure 6 ) show the greatest changes in S0 4 - 
over the region of emission change, with only slight changes 
globally. The individual CTMs and HTAP ensemble mean 
results show that NO x reductions from EA, EU, and NA 
cause both increases and decreases in SO| - over the reduc- 
tion region, which correspond to changes in oxidants (OH, 
H 2 0 2 , 0 3 ); localized decreases in SO 4 - are due in part, 
to localized decreases in OH (auxiliary material Figure S5) 
and decreases in O3. 1 The distributions of H 2 0 2 are not ana- 
lyzed, as they were not reported in the CTM simulations. 
NMVOC reductions from EU increase SO 4 - over north- 
eastern Africa and decrease SO 4 - over western Europe in 
three of the CTMs and the ensemble mean. All four CTMs 
consistently show that the EA NMVOC reduction decreases 
S0 4 - and the SA NMVOC reduction increases S0 4 - 
regionally, while two CTMs show S0 4 - decreases over 
the eastern U.S. from the NA NMVOC reduction. For CO 
reductions from all four regions, the ensemble mean shows 
localized increases, while the individual CTMs differ in the 


’Auxiliary materials are available in the HTML. doi:10.1029/ 
201 1JD017134. 


sign of regional SO 4 - change for all CO reductions except 
SA. Regional increases in SO 4 - from the NMVOC and CO 
reductions can be explained partly by localized increases in 
OH (Figure S5); however, decreases in S0 4 - may be related 
to the differing effects on oxidants (including H 2 O 2 ) and on 
the S0 2 oxidation pathways in each CTM, as discussed in 
the previous paragraph. 

4. Radiative Forcing due to Precursor 
Emission Changes 

[25] Figure 7 (and Table SI in auxiliary material Text SI) 
show the global annual net RF due to 0 3 , CH 4 , and SO 4 - , 
estimated from RTM simulations, first for multimodel mean 
0 3 and CH 4 , and second for multimodel mean 0 3 , CH 4 , 
and SO 4 - , for each SR scenario relative to the base case. 
We calculated the net RF distributions for each perturbation 
scenario (SR2 through SR 6 ) by subtracting (in each grid cell 
and month) the simulated radiative fluxes for the base case 
(SRI) from those for each perturbation. To estimate the 
contribution of the multimodel mean S0 4 - to the global 
annual net RF, we subtracted the net RF results due to O 3 
and CH 4 from the net RF due to O 3 , CH 4 , and SO 4 - for each 
SR scenario, assuming the effects of 0 3 , CH 4 , and SO 4 - are 
additive. To distinguish the contributions of O 3 and CH 4 to 
the global annual net RF, we estimated the net RF due to the 
multimodel mean CH 4 for each SR scenario, using the for- 
mula of Ramaswamy et al. [2001], attributing the difference 
to 0 3 RF. 
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Figure 7. (a) Global annual average RF (mW m 2 ) for the HTAP ensembles of 1 1 models (for O 3 and 
CH 4 forcing) and four models (for SO 4 - forcing) due to multimodel mean changes in steady state O 3 , 
CH 4 , and S 04 _ . Vertical black bars represent the uncertainty in net RF across models, calculated as the 
net RF of the multimodel mean ± 1 standard deviation O3 and CH 4 , for each perturbation (—20% CH 4 
burden, and —20% regional emissions of NO x , NMVOC, CO, and combined), relative to the base 
simulation. The uncertainty estimates for —20% CH 4 account only for the variability in simulated O3 
changes across the CTMs, since all CTMs uniformly reduced CH 4 (1760 ppbv to 1408 ppbv). Vertical 
green bars represent the upper uncertainty bound of SO4 RF across models, calculated as the net RF 
of the multimodel mean +1 standard deviation SO 4 - . The RF of changes in CO 2 uptake by the biosphere 
(yellow) are shown as a range from high to low sensitivity of vegetation to O3, estimated for a single CTM 
(STOCHEM) by Collins et al. [2010]; these estimates are not included in the net RF (supporting data 
provided in Table SI in auxiliary material Text SI). Note the difference in scale between the —20% 
regional (NO x , NMVOC, CO, combined) and —20% CH 4 reduction scenarios, (b) Global, annual average 
RF per multimodel mean change in emissions (mW m 2 /Tg emissions per year) due to multimodel 
mean changes in steady state O3, CH 4 , and SO;* - , and uncertainty (vertical black bars) as the net RF of 
the multimodel mean ± 1 standard deviation changes in O3 and CH 4 per unit change in emissions for each 
perturbation, relative to the base simulation. Vertical green bars represent the upper uncertainty bound of 
SOih RF per unit change in emissions across models. 
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Figure 8. Annual average net RF distributions (mW m 2 ), calculated as the annual shortwave radiation 
minus the annual longwave radiation, due to tropospheric O3, CH 4 , and S0 4 _ for the multimodel mean, 
for each of the precursor reduction simulations (—20% CH 4 burden and —20% regional emissions of 
NO x , NMVOC, CO, and combined) minus the base simulation. Note the difference in scale between 
the —20% regional (NO x , NMVOC, CO, combined) and —20% CH 4 reduction scenarios. 


[26] Computational limitations prevented us from simu- 
lating the RFs individually for each CTM’s 18 SR scenarios. 
Instead, we simulate RFs for multimodel means, and for the 
multimodel mean ± 1 standard deviation O3 and CH 4 and 
the multimodel mean +1 standard deviation S0 4 ~ to account 
for uncertainty in the net RF due to differences in the CTMs. 
We show uncertainty (mean ± 1 standard deviation) for the 
resulting net RF, which includes the uncertainty in O3 and 
CH 4 RFs, as changes in O3 and CH 4 are not strongly cor- 
related among the 1 1 CTMs for most scenarios (Figure S4). 
For NO x reductions, because the RF due to 0 3 opposes that 
of CH 4 , a broader uncertainty range would have resulted had 
we instead simulated together the multimodel mean +1 
standard deviation O3 and the multimodel mean — 1 standard 
deviation CH 4 (and the reverse) to estimate uncertainty. 

[27] Figure 7 shows that O3 and CH 4 RFs have the same 
sign as the tropospheric composition changes in section 3; 
since S0 4 _ is cooling, S0 4 _ RF has the opposite sign. 
NO x reductions from all four regions produce an overall 
positive net RF due to increases in CH 4 , which outweigh 
the negative net RF due to decreases in 0 3 (Figure 7a). 
Negative global net RFs are produced by CO and NMVOC 


reductions, due to 0 3 and CFI 4 decreases, and also by 
the combined precursor reductions, as increases in CFI 4 
from NO x reductions roughly cancel CH 4 decreases from 
NMVOC and CO reductions [Fiore et al., 2009]. The net 
RF due to the combined precursor reduction is 98% to 1 17% 
of the sum of the net RFs (of 0 3 and CH 4 ) due to reductions 
of each individual precursor, across the four regions, 
showing approximate additivity for the different precursors. 

[28] Consistent with the S0 4 changes in Figure 5, NO x 
reductions from EU and SA contribute a positive S0 4 ~ RF, 
while EA and NA NO x reductions produce negative S0 4 _ 
RF. The S0 4 - RFs for NMVOC and CO reductions vary in 
magnitude and sign across the four regions, corresponding 
to the disagreement in S0 4 response across the CTMs 
(Figure 5). We do not estimate the contribution of S0 4 ~ to 
net RF for the combined reductions, since many of these 
perturbations included 20% reductions in S0 2 and aerosols, 
making it difficult to isolate the effect of NO x , NMVOC, 
and CO on SOj~ RF. 

[29] Figure 7 also shows an estimate of the RF due to 
the C0 2 equivalent emission resulting from the influence 
of surface O3 on plants’ ability to remove C0 2 from the 
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Figure 9. Radiative forcing efficiency of O 3 for the 16 SR simulations (SR3 through SR 6 ) for the multi- 
model mean, showing the global, annual average O3 net RF (mW m~ 2 ), calculated as the difference 
between the simulated net RF due to O 3 and CH 4 and estimated net RF due to CH 4 [Ramaswamy et al., 
2001], versus the global, annual average steady state changes in tropospheric 0 3 burden (Tg). The SR 
simulations are distinguished by precursor (color) and region (shape). 


atmosphere, from Collins et al. [2010], based on one F1TAP 
CTM (STOCHEM) and not the HTAP ensemble. The C0 2 
responses to pulse changes from Collins et al. [2010] were 
converted to equilibrium responses by integrating over 
100 years. The range of C0 2 RF represents high to low 
sensitivity of vegetation to O3. With the additional consid- 
eration of C0 2 RF, the global annual net RF due to regional 
NO x reductions changes sign to an overall net climate 
cooling (—0.83 to —4.28 mWm“ 2 for all four regions), while 
the negative net RFs for regional NMVOC and CO reduc- 
tions are reinforced by the addition of C0 2 RF (Table SI 
in auxiliary material Text SI). 

[ 30 ] We normalize the global annual net RF estimates 
(Figure 7a) by each region’s reduction in emissions 
(Table 3). The net RF per unit change in NO x and NMVOC 
emissions is greatest for SA reductions out of the four regions 
(Figure 7b), corresponding to the sensitivity findings in 
section 3.1. For CO reductions, all four regions are approxi- 
mately as effective at reducing global net RF per unit CO 
emissions, consistent with CO’s longer atmospheric lifetime. 

[ 31 ] We compare our ensemble mean global annual net RF 
estimates per unit NO x emissions to those of Naik et al. 
[2005], who used a single CTM and analyzed 10% NO x 
reductions. Our estimates are approximately 32% to 63% 
(EU, NA, SA) and 16% (EA) of those reported by Naik et al. 
[2005], but these differences in net RF (CH 4 and O3 RF 
combined) are small in comparison to the magnitudes of 
CH 4 and 0 3 forcing individually. 

[ 32 ] The net RF distributions (Figure 8) correspond to 
the distributions of total column 0 3 (Figure 3) and S0 4 “ 
changes (Figure 6), where regional to hemispheric RF cor- 
responds to O3 changes and more localized RF is dominated 
by S0 4 _ changes, as illustrated by shortwave forcings 
(Figure S2). Changes in CH 4 influence net RF globally, 


since a uniform CH 4 mixing ratio was specified in each 
RTM simulation. For NO x reductions, we find positive net 
RFs in the SH due to CFI 4 and long-term O3 increases 
globally, but in the NH these positive RFs are outweighed 
by the negative RF of 0 3 decreases (Figure 8). For CFI 4 , 
NMVOC, CO, and combined reductions, negative net RFs 
from O3 decreases in the NH overlay negative RFs globally 
due to CH 4 . While Figure 7 presents globally averaged for- 
cings, the regionally inhomogeneous forcings in Figure 8 
are also relevant for regional climate change [Shindell and 
Faluvegi, 2009]. However, regional RF patterns resulting 
from changes in tropospheric loadings do not directly 
translate to regional climate responses [Levy et al., 2008; 
Shindell et al., 2010], 

[ 33 ] In Figure 9, the relationship between tropospheric 0 3 
burden changes and global O3 RF is strongly linear, giving 
a RF efficiency of approximately 3.27 mW m 2 per Tg O 3 
or 35.6 mWm~ 2 DU -1 (1 DU ~ 10.88 Tg 0 3 [Park et al., 
1999]). This efficiency compares well with those estimated 
in previous studies, 34 to 48 mWm -2 DU -1 [Hauglustaine 
and Brasseur, 2001; Wild et al., 2001; Fiore et al., 2002] 
and 23 mWm“ 2 DU -1 (for NO x ) and 43 mWm 2 DU -1 
(for CH 4 and CO+VOCs) [Shindell et al., 2005], Here 
reductions from EU and NA (with exception of NA NO x ) 
fall at or above the average RF efficiency line, suggesting 
lower RF efficiency. EA and SA reductions (except for EA 
NMVOC) have RF efficiency greater than average, as these 
regions have greater influences on the UT where RF is the 
most efficient [Lacis et al., 1990; Wang et al., 1993; Forster 
and Shine, 1997; West et al., 2009a], 

5. Global Warming Potentials 

[ 34 ] Beyond analyzing RF, the climate impacts of O3 
precursor emissions can be compared with each other, and 
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with the emissions of other species using climate metrics 
such as the global warming potential (GWP). Forster et al. 
[2007] suggest that there are serious limitations to the use 
of GWPs for comparing short-lived species. While other 
metrics have been proposed to compare climate effects, such 
as the global temperature potential (GTP) [Shine et al., 2005, 
2007], none are as widely used as the GWP. We choose to 
analyze the GWP here for comparison with earlier studies. 

[35] The basis for the GWP calculation is the integrated 
RF following a pulse emission. In section 4, O3 RF was 
calculated for equilibrium conditions for the sum of the 
short and long-term 0 3 responses (Figure 7, and Table SI 
in auxiliary material Text SI). Here long-term O3 RF is 
calculated by scaling O3 RF from SR2 by the ratio of 
steady state O3 burden change in a particular SR scenario 
(SR minus base) to those of SR2 (SR2 minus base). The 
short-term O3 RF is then the difference between the steady 
state RF (section 4) and long-term RF. Following Collins 
et al. [2010], RF as a function of time is calculated for a 
one-year emissions perturbation, for each SR scenario. The 
short-term RF components (S0 4 ~ and short-term O3) are 
assumed to be constant over the one-year pulse and then 
drop to zero instantaneously; whereas, the long-term com- 
ponents (CH 4 and long-term O3) respond and decay with the 
multimodel mean CH 4 perturbation lifetime (11.65 years). 
For the 20% CH 4 reduction (SR2), an analytical expression 
is used to calculate the impact of a one-year emissions pulse 
and subsequent decay [ Collins et al., 2010]. This CH 4 per- 
turbation is used to scale the S0 4 - and O3 RFs from SR2 
in Figure 7. The formula by Ramaswamy et al. [2001] is 
used to calculate the CH 4 RF. 

[36] GWPh is given by the RF integrated out to a time 
horizon H and normalized by the change in emissions, 
divided by the equivalent for C0 2 . In Figure 10, and 
Table S2 in auxiliary material Text SI, we present GWPs for 
the 20- and 100-year time horizons (GWP 2 o and GWP 10 o). 
The uncertainties are dominated by the variation in CH 4 
response across the CTMs. For NO x emissions, this uncer- 
tainty is sufficiently large that it is not possible to identify 
the sign of the GWPs. The patterns of GWP100 are very similar 
to the normalized forcings in Figure 7b (the patterns would 
be identical for GWP,*,), whereas the GWP 2 o patterns give 
more emphasis to short-lived O3 and S0 4 - than GWP 10 o- 
For NO x emissions, this brings GWP 2 o proportionally closer 
to zero. 

[37] For NO x reductions, the GWP100 estimates are similar 
to those of Forster et al. [2007] and FuglesFedt et al. 
[2010], though substantially smaller than those provided by 
Shindell et al. [2009]. The O3 contribution from Asian NO x 
found by Berntsen et al. [2005] is within the range of the 
HTAP CTMs’ results, but slightly higher than the HTAP 
multimodel mean. We neglect RFs of nitrate aerosols, but 
Bauer et al. [2007] suggest nitrate contributions to NO x 
GWPs on the order of-80 for GWP 20 and -20 for GWP 100 - 
The NMVOC GWPs (Figure 10) are generally smaller than 
those of Collins et al. [2002]; however, Collins et al. [2002] 
covered a range of individual NMVOCs. For CO, the GWPs 
are comparable to those of Derwent et al. [2001], but smaller 
than those of Berntsen et al. [2005], largely due to the lower 
O3 response of the HTAP multimodel mean; the O3 response 
to Asian CO emissions of Berntsen et al. [2005], however, 


is within the range of the HTAP CTMs. The CO GWP 10 o 
estimates are also smaller than the no-aerosol results of 
Shindell et al. [2009], due to a lower CH 4 response. The 0 3 
contribution to GWP100 for CH 4 is smaller in this multi- 
model study (21% of the direct CH 4 contribution) than the 
25% assumed by Forster et al. [2007], mostly because 
changes in 0 3 above the tropopause are neglected. Adding 
in the contribution of O3 in the lower stratosphere (15% of 
O3 RF) [Sovde et al, 201 1], and that for stratospheric water 
vapor (15% of the CH 4 contribution) would give a total 
GWP100 for CH 4 of 24.2 ± 4.2. 

[38] For NO x and NMVOCs, SA emissions have a larger 
impact than emissions from the other regions. This suggests 
that some latitudinal dependence may be appropriate for 
GWPs of O3 precursors. Note that equatorial or SH emission 
changes were not considered in this study, but Fuglestvedt 
et al. [2010] found a dependence on latitude. European 
NO x emissions have a more negative GWP than other 
regions in the northern midlatitudes, as O3 production in this 
NO x -saturated region is lower. 

6. Conclusions 

[39] We quantify the magnitude and distribution of global 
net RF due to changes in 0 3 , CH 4 , and S0 4 - for 20% 
reductions in global CH 4 and regional NO x , NMVOC, CO, 
and combined precursor emissions. We find that the 20% 
NO x reductions produce global, annually averaged positive 
net RFs, as positive CH 4 RFs outweigh negative 0 3 RFs, 
consistent with previous studies [Fuglestvedt et al., 1999; 
Wild et al., 2001; Naik et al., 2005; West et al., 2007]. For 
CH 4 , NMVOC, and CO reductions, 0 3 and CH 4 RFs are 
synergistic, yielding overall negative net RFs, consistent 
with previous global-scale studies [Fiore et al., 2002; West 
et al., 2007; Shindell et al., 2009]. Including the effects of 
O3 on plant growth and the carbon cycle may change the 
sign of net RF for NO x reductions to an overall net climate 
cooling, in contrast to previous results that neglect this 
effect, while reinforcing the negative net RFs due to 
NMVOC and CO reductions, but future research is needed 
to better quantify this effect. 

[40] By normalizing the net RF estimates by changes in 
emissions (for NO x , NMVOC, and CO), we find that RF is 
more sensitive to NO x and NMVOC emission reductions 
from regions closer to the equator (i.e., SA), consistent with 
our findings that changes in O3 burden per change in emis- 
sions (full troposphere and UT 0 3 for NO x , and UT 0 3 for 
NMVOC) are greatest for SA reductions. RF is more uni- 
formly sensitive to CO emission reductions from each of 
the four regions, which agrees with O3 burden changes 
per unit CO being less variable across the four regions. 
The trends in GWP100 across the four regions, for each 
precursor, reflect the normalized net RF results. Compared 
to GWP, 00, the GWP 20 patterns are influenced more by 
short-term O3 and S0 4 ~. The large uncertainties in the NO x 
GWP estimates, mainly from the variation in calculated CH 4 
responses across the CTMs, limit the determination of 
the sign of NO x GWPs. The estimated GWPs for indi- 
vidual regions are from the largest model ensemble that has 
been analyzed to date, and are broadly comparable to pre- 
vious studies. 
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Figure 10. GWPs for time horizons of (a) 20 years and (b) 100 years for the —20% CH 4 burden 
and —20% regional emissions of NO x , NMVOC, and CO scenarios. The four regions’ estimates 
(labeled “All”) represent the GWP due to the sum of the four regions’ responses (to O 3 , CH 4 , S0 4 _ , 
and all three species [Total GWP]). Uncertainty analysis is as in Figure 7, but also includes the uncertainty 
in the CH 4 lifetimes for the base simulation (SRI) (supporting data available in Table S2 in auxiliary 
material Text SI). 


[41] We find that regional RFs correspond to localized 
increases and decreases in S0 4 ~ burden. 0 3 changes are 
most important for RF on the regional to hemispheric 
scales, and CH 4 influences RF globally, dominating the RF 
response in the SH. The estimated contribution of S0 4 _ 
(direct effect only) to net RF is small compared to the RF 
of O3 and CH 4 . Shindell et al. [2009] found with a single 
CTM that for global NO x , CO, and VOC emissions, changes 
in S0 4 - contributed a RF more comparable in magnitude 
to the RFs of O3 and CFI 4 . Our findings contrast with 
those of Shindell et al. [2009] on the importance of S0 4 “ RF 
due to O3 precursors. However, the robustness of these 
results is limited, since there was substantial variability in 


S0 4 - among only four HTAP CTMs. The effects of O 3 
precursors on S 0 4 via oxidants merit further research, 
using newer models that include improved treatment of 
oxidant-aerosol interactions. 

[ 42 ] We account for variability in 0 3 , CH 4 , and S0 4 - 
across the ensemble of CTMs, but our estimates of uncer- 
tainty only include the variability in CTMs, using a single 
RTM for all RF estimates. We therefore understate the 
uncertainty in net RF. In addition, while we capture the most 
important forcings, a more complete analysis of RF could 
include RFs due to changes in nitrate aerosols (likely 
important for NO x reductions), the indirect effect of aero- 
sols, internal mixing of aerosols, changes in stratospheric O3 
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and water vapor, and changes to the carbon cycle via nitro- 
gen deposition. Finally, we estimate RF due to changes 
in the radiative budget of the global atmosphere, but do 
not estimate the full climate response to regional forcings. 
Future research should link global and regional RF to cli- 
mate responses. 

[ 43 ] Our analysis contributes to the understanding of the 
effects of O3 precursors on global and regional RF, and 
provides motivation for evaluating the climate benefits of 
policies addressing tropospheric 0 3 and its precursors. We 
show that among short-lived O3 precursors, NMVOC and 
CO emission reductions most effectively reduce RF. Our 
GWP estimates could form the basis for regionally specific 
GWPs (or comparable metrics), which may allow O3 pre- 
cursors to be included in future climate agreements and 
emissions trading schemes, as well as provide information 
that could be used to estimate the influence of national or 
regional policies. Variability in GWPs among regions for 
NO x and NMVOCs suggests that regionally specific esti- 
mates would be important. For CO, the consistency in RF 
per unit emissions and GWP 10 o across the four regions 
implies that the error in using a uniform GWP for CO may 
be small. However, the RF per unit CO emissions may differ 
in other regions not studied here. Future studies should 
analyze additional source regions, such as near the equator 
and in the SH, as other precursors show greater sensitivity in 
these regions [e.g., Naik et al., 2005], 
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